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I.  INTRODUCTION 


The  stress  analysis  of  an  arbitrarily  loaded  solid  body  can  only  be 
performed  if  the  stress-strain  relations  are  known  for  an  arbitrary  polyaxial 
stress  state.  In  the  case  of  isotropic,  linearly  elastic  materials,  the  stress- 
strain  law  is  completely  defined  b.v  two  constants,  commonly  Young's  modulus 
and  Poisson's  ratio.  In  ATJS  and  related  graphites,  the  situation  is  far  more 
complicated  because  they  are  neither  linear  nor  isotropic.  In  addition,  there- 
are  some  differences  between  compressive  and  tensile  behavior.  As  it  is 
not  feasible  to  determine  experimentally  the  stress -strain  relation  for  ail 
possible  stress  states,  it  is  necessary  to  combine  a limited  number  of  test 
data  with  theoretical  considerations  in  order  to  develop  a general  stress- 
strain  law. 

One  of  the  most  extensive  sets  of  test  data  for  ATJS  graphite  under 

1 - 5 

combined  stresses  was  obtained  by  J.  Jortner.  These  tests  were  per- 
formed at  normal  and  elevated  temperature  on  both  cylindrical  shells  and 
solid-bar  specimens.  Proportional  loading  tests  were  performed  that  in- 
cluded uniaxial  stresses,  biaxial  stresses,  and  a limited  number  of  triaxial 
stress  states.  Also  included  were  tests  involving  nonproportional  loading. 

Among  the  major  conclusions  from  these  tests  were  the  following: 

1.  The  classical  theory  of  plasticity  does  not  apply  to  ATJS  because 
classical  plasticity  assumes  that  inelastic  strain  is  isovolumetric, 
a behavior  not  encountered  in  this  type  of  graphite.  In  addition, 
the  usual  rule  in  incremental  plasticity  theory  that  plastic  strain 
is  normal  to  surfaces  of  constant  plastic  potential  is  not  obeyed. 
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2.  The  strain  ratios  remain  approximately  constant  during 
proportional  loading  in  compression.  In  the  case  of  tension,  the 
same  is  true  except  for  uniaxial  or  nearly  uniaxial  stress  states. 

3.  ATJS  exhibits  "strain  softening,  11  i.e.  , the  strains  tend  to  be 
larger  under  biaxial  tension  than  under  uniaxial  tension. 

4.  The  tensile  stress-strain  curve  is  always  concave  downward, 
i.e.  , the  tangent  modulus  continually  decreases  with  increasing 
stress.  The  compressive  stress- strain  curve  is  concave  down- 
ward at  stresses  equal  in  absolute  magnitude  or  less  than  the 
tensile  fracture  stress.  In  this  stress  range,  the  behavior  is 
quite  similar  to  the  behavior  in  tension  except  for  Poisson  effects. 
At  much  higher  stresses,  the  curve  flattens  out,  and  there  may 
even  be  an  inflection  point. 

Because  classical  plasticity  theory  was  found  to  be  inapplicable,  the 
main  theoretical  guidance  in  developing  a stress-strain  law  comes  from  a 
considei ation  of  two  limiting  cases.  For  sufficiently  small  stresses,  the 
theory  should  reduce  to  the  appropriate  orthotropic  elastic  theory;  and,  at  any 
stress,  the  equations  should  be  invariant  to  a rotation  of  axes  when  the  aniso- 
tropy is  reduced  to  zero.  A heuristic  account  of  the  development  of  such  a 
stress -strain  law  is  contained  in  the  following  section.  Only  room- temperature 
properties  are  considered,  and  attention  is  limited  to  stress  states  in  which 
no  compressive  stress  is  much  in  excess  of  the  absolute  value  of  the  tensile 
fracture  stress. 
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II.  DEVELOPMENT  OF  STRESS-STRAIN  LAW 


In  the  case  of  a transversely  isotropic  elastic  solid,  such  as  ATJS  at 
very  low  stresses,  the  stress -strain  law  is  essentially  linear  and  can  be 
written  in  cylindrical  coordinates  as:^ 


z v'  v' 

€ z SF  ~ E7  °9  ~ ~E!  °r  ’ 

(la) 

V'  , V 
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r E'  z E 0 E 
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Yr0  G r9  E rG 
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In  the  above  equations,  z is  chosen  in  the  cross-grain  direction,  while  r and  0 

define  the  plane  of  symmetry,  in  which  with-grain  properties  apply.  The 

matrix  relating  stress  and  strain  is  symmetric,  a property  which  follows 

from  the  existence  of  an  elastic  potential.  The  values  of  the  constants  for 

3 

ATJS  have  been  listed  as  approximately: 


= Er  = Eq  = ewg 

= 1.8  X 106  psi  , 

(2a) 

“ Ez  = ECG  = 1,2 

x 10^  psi  , 

(2b) 
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(2d) 


G = G = 8 X 10^  psi 

and 

G'  = G = G Q = 6.  5 X ^05  osi  . (2e) 

zr  20 


Thus, 


—r  = 8.  3 X 10“8  (psi)"1  (2£) 

and 

—•  = 5. 6 X 10"8  (psi)" 1 . (2g) 

So  that  these  results  can  be  extended  into  the  nonlinear  range,  uniaxial 
tension  is  first  considered.  Test  results  that  give  longitudinal  and  trans- 
verse strain  as  a function  of  stress  are  presented  in  Figs.  1 and  2 for  vvith- 

5 

grain  and  cross-grain  loading,  respectively.  It  should  be  noted  that  the 
transverse  strain  is  small  and  approximately  proportional  to  the  loading. 
Consequently,  it  is  assumed  that  in  tension  the  Poisson  effect  behaves 
elastically. 

The  longitudinal  strain  is  nonlineariy  related  to  the  stress,  but  the 
elastic  limit  varies  widely  from  specimen  to  specimen  and  is  in  many  cases 
difficult  to  determine  accurately.  In  cases  in  which  the  nonlinearity  seems 
to  start  imperceptibly  at  very  low  stress,  the  data  can  be  fitted  quite  well 
with  the  assumption  that  inelastic  strain  is  proportional  to  the  square  of  the 
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stress.  A compatible  assumption  for  cases  in  which  an  elastic  limit  stress 
exists  is  that  the  inelastic  strain  is  given  by 


€i  = K¥^- 


(3) 


where  <r®  is  tne  elastic  limit  stress  . When  K is  chosen  to  give,  the  correct 
inelastic  strain  at  fracture,  Eq.  (3)  generally  leads  to  good  agreement  over 
the  entire  stress-strain  curve.  It  turns  out  that  K = 0.  15  leads  to  reasonable 
average  values  for  the  elastic  limit  as  well  as  the  observed  inelastic  strains 
at  fracture,  as  can  be  seen  in  Figs.  1 and  2.  In  practice,  then,  the  elastic 
limit  is  determined  by  requiring  Eq.  (3)  to  be  satisfied  at  fracture  and  solving 
it  for  <r\  The  result  is 


<r 


cr  = 


fracture 


l+6.7(€./e  ) 


e'  fracture 


(4) 


where  ( n and  e ^ are  the  elastic  and  inelastic  longitudinal  strains,  respectively. 
Although  (P  found  in  this  manner  will  not  in  general  exactly  coincide  with  the 
highly  variable  observed  elastic  limit,  the  error  in  total  strain  resulting  from 
the  discrepancy  will  be  so  small  as  to  be  negligible  in  most  cases. 

Next,  the  question  is  addressed:  How  should  this  approach  to  biaxial 
tension  be  generalized?  A plausible  expression  that  reduces  to  Eq  (3)  in 
the  uniaxial  limit  is 


(« 


<r 


x 


(5) 
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where 


F = 


0.  15 


1 


(6) 


when  the  bracketed  expression  is  positive,  and  F = 0 when  it  is  negative. 

The  choice  of  n will  determine  the  amount  of  "strain  softening.  " Table  I 
gives  biaxial  stress  and  strain  data  for  ATJS  at  fracture.  Figure  3 shows 
the  fracture  strains  computed  from  the  experimental  stresses  when  n = 2, 

The  agreement  between  calculated  and  measured  strains  is  quite  good, 
probably  as  good  as  can  be  expected  when  the  variations  'n  properties  between 
specimens  are  c isidered.  In  addition  to  this  empirical  reason  for  choosing 
n = 2,  there  is  a theoretical  reason  that  follows. 

Since  the  biaxial  experiments  just  cited  did  not  include  shear,  they  do 
not  provide  a check  on  whether  shear  should  be  included  in  a general  expres- 
sion for  F.  However,  one  would  expect  that,  as  anisotropy  is  reduced  to 
zero,  F should  be  invariant  with  respect  to  a rotation  of  axes.  An  invariant 
expression  that  reduced  to  a sum  of  squares  of  the  direct  stresses  in  the 
principal  axis  system  is 


J 


2 

1 


+ 2J2 


+ 


2 2 

i^j  7ij 


(7) 


where  Jj  and  are  the  first  and  second  stress  invariants,  respectively. 
Thus,  a plausible  expression  for  F in  the  general  case  becomes 
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Table  I.  Jortner's  Biaxial  Failure  Data 


<r  /<ra 

z 0 

cr  (ksi) 
z 

<r0  (ksi) 

c 

z 

€e 

1:0 

3.940 

0 

0.00548 

-0.00038 

1:1.26 

3.  358 

4.  246 

0.00462 

0.00348 

0:1 

0 

4.  786 

-0.00033 

0.0039 

-1:1 

-4. 120 

4.220 

-0.0063 

0.0043 

£ 

Each  value  listed  is  the  average  of  three  to  five  experiments. 
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Fig.  3.  Comparison  of  Measured  and  Computed  Fracture 
Strains  Under  Biaxial  Loading 
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when  the  bracketed  expression  is  positive,  and  F = 0 when  it  is  negative. 

With  this  definition  of  F,  the  general  stress- strain  law  for  tension 
becomes: 


o 

:=#-  (1  + F)  - 8.  3 X 10"  (<rQ  + <r  ) 
I:  Or 


,-8  , °0 


J * T—N 

z E 


,-8 


■\4 


(9a) 

10-8<r 

(9b) 

r 

(1  + F)  , 

(9c) 

V0z“  (1  + 


Yrz=  (I  tFI 


rz 


and 


\'ro  s H + F)  G 


(9d) 


(9e) 


(9f) 


Note  that  the  values  for  E and  G are  given  in  Eq.  (2).  The  constants  needed 
to  determine  F were  found,  in  fitting  the  tensile  curves  of  Figs.  1 and  2,  to 
be 


I 


ft 
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(10a) 


<r  =0.82  ksi 
z 


and 


o-  0 = c-.0  = 1.22  ksi 
r 0 


(10b) 


It  was  found  that  the  compressive  stress- strain  curve  was  close  to  the 
tensile  curve,  but  that  the  Poisson  effect  was  quite  different.  In  compression, 
the  transverse  strain  was  proportional  to  the  longitudinal  strain,  rather  than 
to  the  longitudinal  stress.  Accordingly,  under  combined  compressive 
stresses,  the  stress-strain  law  is  assumed  to  be: 


€ = (1  + F) 
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% s 
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The  only  difference  between  the  equations  for  tension  and  those  for 
compression  lies  in  the  off-diagonal  terms  that  give  the  Poisson  effect.  In 
mixed  states  involving  both  compression  and  tension,  the  Poisson  term 
associated  with  a compressive  stress  will  be  calculated  as  in  Eq,  (11),  while 
the  Poisson  term  associated  with  a tensile  stress  will  be  calculated  as  in 
Eq.  (9). 


III.  COMPARISON  WITH  TRIAXIAL  STRESS  DATA 


The  stress-strain  laws  just  formulated  were  based  on  uniaxial  tests  in 
compression  and  tension,  and  on  general  theoretical  considerations.  Their 
form  implies  superposition  of  Poisson  effects.  While  such  superposition  is 
a requirement  when  Hie  relations  are  linear,  it  cannot  be  assumed  to  be  valid 
in  the  nonlinear  case  with  which  this  paper  is  concerned.  It  therefore  can 
only  be  used  with  confidence  if  experimentally  confirmed.  Such  confirmation 
was  shown  in  the  preceding  section  in  the  case  of,  biaxial  fracture  data  in 
tension.  Now  the  equations  will  be  checked  against  triaxial  proportional 
loading  experiments  in  which  the  stresses  in  the  plane  of  symmetry  are  equal 
compressions  and  the  stress  in  the  cross-grain  direction  can  be  either  tensile 
or  compressive  in  nature. 

The  test  data  from  the  triaxial  tests  are  shown  in  Fig.  4,  which  Is  taken 
from  the  reference  cited  in  footnote  5.  Superposed  on  the  experimental  curves 
are  calculated  points  based  on  the  use  of  Eqs.  (9)  and  (11).  The  good  agreement 
between  theory  and  experiment  appears  to  validate  both  the  above-noted  neces- 
sity of  two  different  treatments  of  the  Poisson  effect,  depending  on  whether 
the  associated  stress  is  positive  or  negative,  and  the  assumption  of  superposi- 
tion of  Foisson  effects.  It  also  appears  to  indicate  that  the  nonlinear  coupling 
between  stress  components,  which  is  assumed  in  Eq.  (8),  is  at  least  approxi- 
mately correct.  In  making  this  comparison,  it  was  noted  that  the  specimens 
were  bars  and  that  bars,  were  slightly  stiffer  in  the  cross-grain  direction  than 
the  cylinders  to  which  the  uniaxial  laws  were  tailored.  To  account  for  this, 

«•  ® was  taken  to  be  1 . 04  ksi . 
z 
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Fig.  4 
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. Across-Grain  Stress-Strain  Responses,  Triaxial  Tests, 
Compared  with  Theory  (adapted  from  reference  cited 
in  footnote  5) 
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IV.  NONPROPORTIONAL  LOADING 


All  the  data  discusseo  up  to  this  point  were  obtained  from  proportional 
loading  tests,  and  the  equations  developed  to  represent  them  are  in  finite 
rather  than  incremental  form.  Such  a form  is  appropriate  for  proportional 
loading,  but  if  the  inelastic  strains  associated  with  a given  stress  state  depend 
significantly  on  the  loading  path,  the  equations  just  developed  cannot  be  used 
for  nonproportional  loading. 

Nonproportional  loading  paths  and  the  corresponding  strains  are  shown 
in  Figs.  5 and  6,  respectively,  taken  from  the  reference  cited  in  footnote  5. 

The  fact  that  three  different  loading  paths  to  the  same  final  stress  state  led  to 
three  different  strain  states  suggests  some  path  dependence,  but  the  evidence 
is  really  not  clear.  Incremental  plasticity  theory  would  predict  that,  when  a 
specimen  is  first  loaded  in  one  direction  and  then  a load  is  added  at  right 
angles,  the  strain  in  the  first  direction  will  be  bigger  than  if  the  loads  are 
applied  in  the  reverse  order.  This  qualitative  prediction  is  fulfilled  in  the 
case  of  c but  not  in  the  case  of  < Straightforward  application  of  the 
equations  of  the  present  paper  to  path  bx  of  Fig.  5 leads  to  the  data  represented 
by  the  triangles  shown  in  Fig.  6.  It  is  evident  from  a comparison  of  theory 
and  experiment  for  point  b that  specimen  87  was  somewhat  weaker  in  the 
cross-grain  direction  than  assumed  in  the  theory,  which  could  account  for  the 
discrepancies  in  the  remaining  points.  In  any  event,  it  is  possible  to  show 
from  very  general  considerations  (omitted  here  for  the  sake  of  brevity)  that 
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Fig.  6.  Strain  Trajectories  (as  given  in  reference  cited  in 
footnote  5)  Compared  with  Theoretical  Points 
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incremental  theory  would  lead  to  larger  discrepancies  with  the  experiment 
them  would  application  of  the  present  theory  without  consideration  of  the 
loading  path.  The  latter  procedure  therefore  appears  to  be  indicated  until  a 
more  complete  understanding  of  the  path  dependence  question  is  available. 


V.  RELATION  TO  OTHER  STRESS-STRAIN  LAWS 


Three  polyaxial  stress-strain  laws  have  previously  been  proposed  for 
transversely  isotropic  graphites.  Here  they  will  only  be  identified  and  briefly 
discussed,  since  a detailed  comparison  with  the  current  approach  and  with 
test  data  is  beyond  the  scope  of  this  paper. 
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The  first  stress-strain  law  is  due  to  Weng.  It  states,  in  effect,  that 
all  strain  components  increase  proportionally  during  proportional  loading,  an 
assumption  that  appears  to  be  valid  for  compression,  but  only  approximately 
true  for  tension.  His  function  p,  which  corresponds  to  our  (1  + F),  was  not 
formulated  analytically  and  did  not  include  contributions  from  shear  stresses. 
His  results  therefore  do  not  constitute  a general  stress- strain  law,  and  in 
addition  they  seem  not  to  be  in  as  good  agreement  with  the  experiment  as  the 
present  formulation. 

An  incremental-type  theory  patterned  after  classical  plasticity  theory 

8 

was  proposed  by  Weiler.  The  theory  explicitly  treats  thermal  as  well  as 
load-induced  stresses  and  assumes  a bilinear  or  trilinear  stress- strain 
relation.  However,  the  present  author  believes  that  the  experiments  cited  in 
footnote  5 invalidate  the  incremental  plasticity  theory  for  ATJS-type  graphite. 
The  use  of  incremental  plasticity  theory  is  therefore  unwarranted,  particu- 
larly when  the  additional  complexity  of  incremental  theories  is  taken  into 
consideration. 

g 

Jones  and  Nelson  have  proposed  a deformation-type  theory  in  which 
the  various  material  parameters  appearing  in  Eq.  (1)  are  independently 
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variable  with  the  stress  and  strain  level.  The  rotational  invariance 
requirement  is  met  by  relating  the  values  of  those  parameters  to  the  contracted 
matrix  product  of  the  stress  and  strain  tensors.  This  approach  appears  to  be 
a good  deal  more  complicated  to  apply  than  that  advanced  in  this  paper.  In 
the  author's  view,  the  success  of  the  present  approach  makes  the  additional 
complications  questionable,  especially  as  it  is  not  clear  that  any  greater 
accuracy  would  result. 
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VI.  DISCUSSION 


This  paper  presents  a new.  nonlinear,  polyaxial  stress-strain  relation 
for  ATJS-type  graphites  that  r educes  to  the  known  elastic  relations  at  low 
stresses  and  is  in  good  agreement  with  extensive  biaxial  and  triaxial  data  at 
higher  stresses.  It  is  a deformation-type  theory,  which  means  it  is  strictly 
applicable  only  to  proportional  loading.  However,  it  gives  a rather  good 
approximation  (better  than  simple  classical  incremental  theories)  even  in 
nonproportional  loading.  The  author  believes  that,  although  this  theory  was 
formulated  to  apply  specifically  to  ATJS,  with  minor  changes  in  moduli  and 
Poisson's  ratios  it  can  also  be  used  for  related  graphites. 

As  written,  the  equations  give  strains  as  functions  of  the  stresses. 
Because  of  the  nonlinear  coupling  between  the  stresses,  an  exact  analytical 
formulation  of  the  stresses  in  te-ms  of  strains  is  not  practicable.  However, 
the  problem  can  readily  be  solved  iteratively  as  follows.  Given  the  strains, 
an  arbitrary  numerical  value  for  F (such  as  0.  5)  can  be  assumed.  This 
converts  the  equations  into  a set  of  simultaneous  linear  algebraic  equations 
that  can  readily  be  solved  for  the  stresses.  With  these  stress  values,  a 
corrected  value  for  F can  be  computed  and  the  process  repeated,  continuing 
until  a satisfactory  degree  01  convergence  is  achieved. 

In  summary,  the  analytical  formula  for  the  polyaxial  stress-strain 
relations  of  ATJS  proposed  here  is  a generalization,  guided  by  theoretical 
considerations,  of  a large  amount  of  test  data.  The  result  is  neither  rigorous 
nor  unique,  but  it  is  simple,  surprisingly  accurate,  and  relatively  easy  to  use. 
It  should,  therefore,  be  of  value  in  structural  analysis. 
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LABORATORY  OPERATIONS 


The  Laboratory  Operation*  of  The  Aerospace  Corporation  ia  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opment* is  vital  to  tie  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-powcr  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmonpheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optica;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory;  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Physics  Laboratory;  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  rnagnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  system*. 
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